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Abstract 

Recently it was argued that gravity with the squire of the Ricci tensor can be stabilized by adding 
constraints to the theory. This was so far demonstrated for fluctuations on the Minkowski/de Sitter 
background. We show that the same scheme works equally well for removing Ostrogradski’s ghost 
from fluctuations on a cosmological background in generic f{R, R'fiy, C'^j^p(^)-type theories of gravity. 
We also derive the general formula for the spectrum of primordial tensor perturbations from the 
stabilized theory. 
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I. INTRODUCTION 


Physical laws describing the time evolution are written in the form of differential equations 
up to second order in time. This is indeed the case for classical mechanics, Maxwell’s theory 
of electromagnetism, and general relativity. Then, what if one has the evolution equations 
of higher order in time? The answer is tragic; one will encounter Ostrogradski’s ghost that 
is the generic instability in non-degenerate higher derivative theories [1-3]. Ostrogradski’s 
instability can be illustrated by the following simple example [4]: 


S = 



- Vil) 


( 1 ) 


This action yields the fourth-order equation of motion for q. Dehning the canonical coordi¬ 
nates Qi := q, Q 2 '■= q and their conjugate momenta Pi, P 2 , one obtains the Hamiltonian, 

H = PiQ2 + ^ + V{Qi). ( 2 ) 


This Hamiltonian linearly depends on Pi and hence is not bounded from below, signaling 
the instability. 

General relativity, a second-order theory for the metric, is a healthy theory from this 
viewpoint. Nevertheless, it is sufficiently reasonable to consider gravitational theories be¬ 
yond general relativity. Higher powers of the curvature tensors such as and 
are expected to be the low-energy manifestation of quantum gravity. In recent years phe¬ 
nomenological modification of general relativity has been studied extensively in order to 
account for the present accelerated expansion of the Universe, which involves an arbitrary 
function of the curvature tensors or a non-minimal coupling to the scalar held. Due to the 
higher derivative nature, those theories are often plagued by Ostrogradski’s instability. 

There are several ways to evade this instability issue. From an effective held theory point 
of view, if one takes the UV cutoh below the scale at which the ghost emerges then the 
instability is not necessarily problematic and thus could be circumvented. More directly, 
one would consider the class of theories that violates the assumptions of Ostrogradski’s 
theorem. The R? model [5] and its /(P) generalization [6, 7] are degenerate, and hence are 
free of ghosts. The Galileons [8-10] and the Horndeski family of scalar-tensor theories [11- 
13] have manifestly second-order held equations despite the higher derivative nature of the 
Lagrangian, and therefore we do not need to care about Ostrogradski’s ghost. (In fact, it is 
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well known that f{R) theories can be recast in a scalar-tensor theory that is in a subclass 
of the Horndeski theory.) Recently, an approach in this direction has been pushed forward 
and healthy theories beyond Horndeski have been developed [14-20]. 

Yet another way of getting around the ghost is adding constraints to the theory to reduce 
the dimensionality of the phase space, as proposed in Refs. [4, 21]. In Ref. [21], the theory 
described by [22] 

A/p f 

S = ^ (R - 2A + aR^ + (3) 

was examined. Although this theory contains unstable modes as it is, it can be made stable 
by the appropriate addition of constraints at least at the level of the quadratic action for 
fluctuations. This was demonstrated for the fluctuations on the Minkowski and de Sitter 
backgrounds [21]. 

The purpose of the present paper is to extend the work of [21]. We consider the more 
general action of the form 

S = j dSv^/(R, (4) 

where / is an arbitrary function, and stabilize the quadratic action for the fluctuations on 
the Friedmann-Lmaitre-Robertson-Walker (FLRW) background. Note that since the Weyl 
tensor C^upa can be expressed in terms of the Riemann tensor, the Ricci tensor, and the Ricci 
scalar, the above theory is nothing but f {R, gravity. This class of gravitational 

theories has been explored extensively [23-30]. As an application, we derive the general 
formula for the primordial power spectrum of the tensor modes from the stabilized action. 

The paper is organized as follows. In the next section we give a quick review on the insta¬ 
bilities in f{R, C‘f^ypa) gravity. We then stabilize the quadratic action for cosmological 
perturbations by adding constraints in Sec. III. In Sec. IV, we evaluate the amplitude of 
primordial tensor modes in the stabilized theory. Section V is devoted to discussion and 
conclusions. 

II. INSTABILITIES IN /(R, R^^, GRAVITY 

We start with reviewing how ghost instabilities appear in the theory described by 
the action (4). The background we consider is given by the FLRW metric, ds^ = 


3 



a^{ri) (—+ 6ijdx^dx^). Let us use the notation 


r,:=R, := (5) 

and write /o := /, /* := df/du, fij := d'^f/dridrj, and fijk := f/drtdrjdrk evaluated 
at the background, he., ri = 6 {R} + R') /a^, r2 = 12 + R^R' + /a^^, and ra = 0, 

where R := a'/a and the prime denotes differentiation with respect to the conformal time 

7 ]. 

The background equation can be derived simply by substituting the metric ds^ = 
—N'^{ri)dri‘^ + a^{ri)6ijdx’‘dx^ to the action (4), varying it with respect to N and a, and 
then setting N = a. Variation with respect to N gives 


0 = aVo - da^R'h - 12 [( 2/2 + 3/n) R {2R^ - R”) + / 2 ^' {R^ + 2R')] 

144 r o 

-fi2R 3R^ - R^R" + R {R'f + 2R^R' - 2R'R'‘ 


144 


/ 22 'H {R^ + 2R') AR^ - R^R” + 2R {R'f - 2R'R 


( 6 ) 


This equation will be used in the following calculations to eliminate or rearrange some 
coefficients in the quadratic actions for fluctuations. Variation with respect to a gives another 
cumbersome equation which is not used in the present paper. 


A. Tensor perturbations 


Let us first look at Ostrogradski’s instability of tensor perturbations. It is convenient to 
parametrize the tensor perturbations as 


ds^ = 


-di] + ( 5ij + hij + -hikhkj ) dx'dx^ 


(7) 


as this dehnition yields at quadratic order in fluctuations. Substituting this 

metric to the action (4), we obtain the quadratic action for the tensor perturbations: 
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This can be simplified by integration by parts as 



where we defined 


{Ahfj + Chijd'^hij) + ^ h”j + 2h[jd‘^h[j + {d'^hij) 


(9) 


/ 3 := 2/2 + 4 / 3 , 


( 10 ) 


and 

A := 2/i + ^^{n^ + W) - 4/'^, (11) 

C := 2/i + + H') + - 4 (/2 - 2/^'). (12) 

a"' ' 

Throughout the paper we assume that /3 = P{ri) never vanishes.^ The resultant action 
turns out to be of the same form as that on the Minkowski/de Sitter background in the 
aR^ + l3Rfj,yR^’^ theory [21], though in the present case the coefficients A, C, and [3 are 
time-dependent in general. The quadratic action (9) contains higher derivative terms, and 
hence Ostrogradski’s instability appears as expected. 

We confirm the presence of Ostrogradski’s instability in a more rigorous way by the 
Hamiltonian analysis. Defining the canonical coordinates as 

fty = Ay ^ 7r« = i [oMAA - n (A'" - 28%) - I3’h",] . (13) 

?y = AA ^ p« = |a" , (14) 


we obtain the Hamiltonian 


H = d^x 
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jAy (d^CS^ + pcdcd) Ay - by (o^A + 2dS^) % 


. (16) 


This Hamiltonian linearly depends on TTjj as can be seen in the first term, implying that the 
Hamiltonian is not bounded from below. Thus, we see that the higher derivative nature of 
the action (9) gives rise to the ghost instability. 


^ If / depends on the curvature invariants through the Gauss-Bonnet combination, Q := —4:R^^ +, 

i.e., f = f{R,G), we identically have /3 = 0. Therefore, we do not consider the f{R,G) class of theories 
in this paper. Linear cosmological perturbations are healthy in f{R,Q) gravity [31], though it was found 
in the end that ghost degrees of freedom cannot be avoided on less symmetric backgrounds [32]. 
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B. Vector perturbations 


The vector sector of the metric perturbations is given by 


ds^ = [— (l — BiB^) d?7^ + 2Bidx^dr] + 5ijdx^dx^^^ , (16) 


with diBi = 0. Here we added in the (00) component so that = a^. Substituting 
the metric to the action (4), we obtain the quadratic action for the vector perturbations: 



where f3 and A are the same as the corresponding quantities defined in the previous sub¬ 
section, and Vij := diBj. We thus see that in the theory (4) the vector perturbations are 
dynamical in general. However, the vector sector is free of any instabilities provided that 
/9 > 0 and A < 0. Note that here again the quadratic action takes the same form as that 
on the Minkowski/de Sitter background in the aR^ + fdR^yR^'^ theory [21], but with the 
time-dependent coefficients. 

In terms of the canonical momentum conjugate = {f3/2)v^^, the Hamiltonian of the 
vector sector is written as 


H = d^x 




a^A 




(18) 
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This Hamiltonian is bounded from below and therefore the vector sector is stable for (3 > d 
and H < 0. 


C. Scalar perturbations 

To simplify the manipulation we fix the gauge inside the action. It is probably the most 
suitable to take the flat gauge, and it is indeed possible to do so at the action level [33]. The 
metric can thus be written as 

ds^ = a? [—(1 -I- 2^)drf‘ + 2diBdr]dx^ + 5ijdx^dx^]^ . (19) 

The quadratic action for the scalar perturbations is given by 

= y d^x [bo ($0^ + hi + B'Y + + ho^B' 

Tcid)^ + C2<ha^<h + c3<hB + 04^2 + , (20) 
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where B := d^B. The coefficients in the above action are 


ho 

:= , 

(21) 

hi 

;= -X + - , 

3 3 

(22) 

h2 

;= AHX , 

(23) 

h. 

:= A{2H' -H^)I , 

(24) 


where 

X := 2/2 + 3/n + + ‘2n') /12 + + 2nf f22 • (25) 

a'* 

The concrete expressions for q {i =1-5) are lengthy and are summarized in Appendix A. 

It is not obvious from (20) whether or not the scalar sector contains unstable degrees of 
freedom. Let us therefore take a careful look at the Hamiltonian for the scalar perturbations. 
We choose to use (<h, B) as canonical coordinates, and then the corresponding canonical 
momenta are given respectively by 


The Hamiltonian is 
H = f d^x 


7r<j) — 26o*h^ + ^2^^, 

TTs = 2bi (^2$ + B') + b2^' + 53 $. 




4/S 




-Ci<f)2 - C2$a2<h - C3$H - C^B^ - c^Bd^^ 


To see the instabilities it is more convenient to perform a canonical transformation, 


7r<j> 

TTs 

$ 

B 


7r$, 

TTg - 




$ -I- 

sn' 


= B. 


The generating function F of this transformation is given by 

F= t i^x + ^) + JbB - X + 2hX&‘B) 


(26) 

(27) 


(28) 


(29) 

(30) 

(31) 

(32) 


(33) 
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The above canonical transformation turns the Hamiltonian into 


H = 



TTg — $63 — 26i(9^$) + 


TTS 


2X + /3 




aM <h - 


B 

m 


+ •■■ ,(34) 


where we have written only the terms that are relevant at high momenta. This Hamiltonian 
clearly shows how instabilities arise: it is required for stable kinetic terms that /3 > 0 and 
X > 0, but then the gradient instability occurs at high momenta, as seen in the coefficient 
of the term. 


III. STABILIZATION OF /(ii, , (72 ) GRAVITY 

In the previous section we have reviewed instabilities of cosmological perturbations in 
/(/?, gravity. For all types of perturbations, it is found that the basic structure 

of the quadratic actions is very similar to that on the Minkowski/de Sitter background in 
aB^ + gravity [21]. This suggests that one can stabilize the theory, in the same 

way as in Ref. [21], by imposing suitable constraints while maintaining the renormalization 
properties. In this section, we demonstrate that this is indeed true. 


A. Tensor perturbations 


Following Ref. [21], we introduce a auxiliary tensor held into the action (9): 
S = {Ahf. + Ch,jd‘^h,,) 


+ 






2 


Variation with respect to Xij gives the constraint 


(35) 


Xij — h'X + 2d^ hij — 0 


(36) 


Substituting this constraint back into the original action (35), we arrive at 

S=^jd^x {h[j {a^A - 2(3d^) hd + hij [{a^C + 2^") - 3/3d^d^] h^j} . (37) 

The dangerous second time derivative h'X can thus be removed from the action while retaining 
higher spatial derivatives acting on h^j and hX. (The spirit here is the same as that of Hofava 








gravity which is power-counting renormalizable [34].) The reduced action implies that the 
tensor sector becomes free of Ostrogradski’s instability. 

To conhrm the stability let us construct the Hamiltonian of the tensor sector. The 
canonical coordinates we choose are 


= fty M J (ft” - Aft - 2d%) - f)'(ft" - A.^)] , 

^ ftft ^ p« = I (ft" - Ay), 

Ay = Ay < -^ y" ^ 0. 


(38) 

(39) 

(40) 


The Hamiltonian is then given by 


H = 




TT 


V 


+ ^Pij + (PiJ - 2^ K 


/S 


/S 


^hij {a^Cd^ + jSd'^d'^) hij - ^qij [a^A + 2/38'^) qij 


(41) 


The primary constraint is =0. The consistency of the constraints generates the following 
set of secondary constraints: 


Pij 


Ci2 


8 hij ^ 0 , 


1 2 A 

TTij - -a Aqij ^ 0, 


[a^A — 2(38'^) Xij —a^Apij — (^a^C8‘^ {38'^8'^) hij ^ 0, 


(42) 

(43) 

(44) 


where ~ stands for weak equality. These are the second class constraints. We use these 
constraints to eliminate {Xij^'wD and {qij,pij), and obtain the reduced Hamiltonian, 

Hr = I d^x -^^Tiii (a^A — 208'^) TTii + Xhii (—a^C8'^ + 308‘^8‘^) hij , (45) 


which is positive definite if 


0{p) > 0, A{p) > 0, C{p) > 0. 


(46) 


This is the sufficient conditions for the absence of instabilities. Note, however, that if C 
becomes negative for a sufficiently short period then the tensor sector is still stable because 
only low momentum modes develop instabilities whose time scales are bounded from below. 

Since each tensor variable has two independent components, the original theory (9) con¬ 
tains eight degrees of freedom in phase space. We then add four second class constraints, 
leaving four degrees of freedom in phase space in the constrained theory. The ghost modes 
in the tensor sector can thus be removed. 
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B. Vector perturbations 


The vector perturbations are not harmful in itself; as we have seen in the previous section, 
the vector modes are stable for /? > 0 and A < 0. However, this is incompatible with the 
stability conditions for the tensor perturbations (46). For this reason, we are going to remove 
the vector modes from the theory. 

Introducing a auxiliary vector field Aj, we consider the modified quadratic action 

a^A , 


^ = i I ^ 


^ij) A Vijd Vij 


+ 


—f > 


with Xij 


:= diXj. The canonical coordinates are chosen to be 


/S 


Vij = Vij < -^ TT*^ = - () , 

Xij = Xij <—)■ d!;' = 0, 


(47) 


(48) 

(49) 


and the Hamiltonian is 


H = / d^x ( ^ + nijXij - jVijdHij - 


We can derive secondary constraints from the primary constraint tt? = 0 as 


(50) 


TT*-^ 0, 


[a^A + (38“^) Vij ^ 0 , 

(oM + pd'‘) (a,, + 


(51) 

(52) 

(53) 


Substituting these constraints back to the Hamiltonian, we see that the reduced Hamiltonian 
vanishes, indicating that there are no vector degrees of freedom. 

C. Scalar perturbations 

Along the same line as the stabilization procedure on the de Sitter background [21], we 
introduce a auxiliary scalar field A to modify the quadratic action as 

S = J d8x[boid>'f + bi + B' - Xf + b2^'{B' - A) + b^^{B' - X) 

+Ci$^ + C24>a2$ + c3<hH + cJ3‘^ + . (54) 


10 








The canonical momenta are now given by 


7r$ = 26o$' + fc2(S'-A), 

ITS = 26, + B' - A) + 62 ®' + 63 *, 

tta ^ 0 , 


(55) 

(56) 

(57) 


and the Hamiltonian is 

H = f d^x 




-Ci$2 - C2<h<92$ - c3<hH - c4H^ - c5Bd^d) 


(58) 


The primary constraint reads tt^ = 0, and the consistency of the constraints generates the 
following secondary ones: 


ttb ~ 0 , 

C 3 $ + 2 C 4 H + C5(9^$ 0 , 

C 4 A T ■ ■ ■ ~ 0, 


(59) 

(60) 
(61) 


where the last equation fixes A. Substituting these constraints to the Hamiltonian (58), we 
obtain 


Hr = I d^xl (71^ + 


V ^ 2I + /3 

C 3 C 5 


C3 363 

4(2I + /3)_ 




- C2 - 


63 $ + —(r$)^ 

2C4 J 4c4 ^ 


(62) 


We thus see that only the two degrees of freedom are left in phase space. Performing a 
canonical transfromation, 

QHI 


TT^, — TT^, 


$ = $, 


21 + ^ 


( 63 $ + 261 ^ 2 $) , 


whose generating function is given by 

F = J d^x 

we obtain 


7r$4> - + 26i4><92<h) 


2 X + /3^2 ^ ^~n 2 

24-H2j^^-3> + 


Hr = I d^x + di<^^ - d2^d‘^^ 


(63) 

(64) 

(65) 

( 66 ) 
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where 


hHxy 

21 +/ 3 ; 


c| 3feo 

" 4^: “ + 4(21 + ff) 


-3 


(^2 — C2 


C3C5 

2C4 


&3 + 2CHX)'. 


(67) 

( 68 ) 


The stability of the tensor sector has already imposed /3 > 0. For a stable kinetic term we 
therefore require that 


X > 0 or 2X + /3 < 0. (69) 

Requiring that C 4 > 0, di > 0, and (^2 > 0 is sufficient for Hr to be positive dehnite. 
However, one may relax the condition and allow for negative di and ^ 2 , as what is crucial 
is the time scale of instability growth. In light of this, we have to avoid encountering the 
rapid, catastrophic growth of the gradient instability at high momenta, so that at least we 
must require C 4 > 0 , while di and d 2 can be negative for a sufficiently short period. 


IV. PRIMORDIAL TENSOR SPECTRUM IN HIGHER DERIVATIVE GRAVITY 
WITH CONSTRAINTS 

The quadratic action for the tensor perturbations with constraints (37) has a non¬ 
standard kinetic term as well as a higher spatial derivative term. It would therefore be 
interesting to explore whether or not this novel structure of the quadratic action gives rise 
to characteristic imprints on the primordial tensor spectrum from inflation. To do so we 
work in a (quasi-)de Sitter background without assuming any particular form of the function 
During inflation we have "H/a ~ const and H' ~ H?, which leads to /3 ~ 

const and 

A ~ C ~ 2/1 -I- 16 / 2^7 ~ const. (70) 

Thus, ignoring slow-roll suppressed contributions, it suffices to consider 

S = ^ J d^x [h[j (a^A - 2f3d‘^) h[ - + hij {a^Ad'^ - SPd'^d'^) hij] , (71) 

with a = {—Hri)~^, where H is the Hubble scale during inflation. 

Following the standard procedure of quantization, we move to the Fourier space and 
introduce the canonically normalized variable 

Vij{rj]k) := ]^[a^A + 2l3k'^y^‘^hij{rj]k). (72) 
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The equation of motion is given by 


f'' + ^ 0 , 


(73) 


where 


with 


2 l + S^k'^r]'^ 

:= 


1 + 2^^2772 


k^- 


1 + 2^k?rf 7f ) 


i--= 


A 


(= const). 


Since ~ (3/2)/c^ as \kr]\ —)■ cxo, the appropriate initial condition is 


7;. . ~ 

Uij — 




(3/2)V4A;1/2 

The equation of motion (73) can be written using y := —kr] as 

2 1 _ n 

“ (i + 2e2/2)7/2j'^*^ “ • 


d Vij 1 + 3^7/ 


(74) 


(75) 


(76) 


(77) 


d7/2 1 + 2^7/2 

The growing solution for 7 / 1 is given by Vij cx I/ 7 /. Therefore, the super-horizon solution 

to Eq. (77) with the initial condition (76) is of the form 

1 C(0 


7; . . 


^ , (78) 

y 

where C(.^) should only be characterized by The primordial tensor spectrum is then given 
by 

27/2 


Vk = 


(79) 


A 7r2 

The power spectrum is nearly scale-invariant, but the amplitude is modified from the stan¬ 
dard result by the model-dependent factor C'^/A. 

We numerically solved Eq. (77) for different values of ^ to fix C(.^). From Fig. 1 it is 
found that 

1 




(80) 


1 + 4’ 

where s is nearly constant even for larger ^ than plotted and s > 2.4. Thus, the primordial 
tensor spectrum is evaluated as 

1 2772 


Vh ^ 


(81) 


A + sjSH^ ir2 

It is worth emphasizing that this result is obtained without assuming any particular form 
of the function /(72, 
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FIG. 1: 1/C^ as a function of 

V. DISCUSSION AND CONCLUSIONS 

In this paper, we have extended the work of Ref. [21] to more general theories of gravity on 
a less symmetric background, and shown that cosmological perturbations in f{R, Rfiuprr) 
gravity can be stabilized by adding constraints to the theory at the level of the quadratic 
action. We have found that the stabilized theory has two tensor and one scalar polarizations. 
This indicates that the propagating degrees of freedom are the same as those in scalar-tensor 
theories, though it involves safe higher spatial derivatives of the metric fluctuations such as 
dtd'^Sgpiy and d’^d'^Sg^u- R would therefore be intriguing if one could identify the correspond¬ 
ing scalar-tensor theory. We expect that the ADM description of scalar-tensor theories in 
the unitary gauge is the optimal way for this purpose, following and generalizing the recently 
developed approach toward single-scalar theories beyond Horndeski [14, 15]. Helpful hints 
for guessing the corresponding scalar-tensor description would be obtained by going beyond 
the quadratic action or by examining perturbations on more general backgrounds such as 
Kasner spacetime. We hope to report our developments in this direction soon. 

We have focused on the universal structure of the action for the tensor perturbations 
during inflation, and derived the power spectrum of primordial gravitational waves from the 
stabilized theory. It would be interesting to evaluate the power spectrum of the curvature 
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perturbation to confront the stabilized theory with observations. We will come back to this 
issue in a future publication. 
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Appendix A: Coefficients a in the quadratic action 

The following is the list of the coefficients q in the action for the scalar perturbations: 
Cl := + 36/ii \aH^ - 2HH" + {Uf + H^H'] + 24/2 - 2HH" + 


72 


fi2 22-H® + 4 {H'f + 1977^77' - 877^77" - imU'H" + 1177" (77 


//^2 


216 


/iii77 (77" + 77') (277^ - 77") 


/22 1 277'' + 2 (77')^ + 2l7fH' - + 1177^ (77 


6-17' 






1077 {H'f 77" + 1177" {H'f - 1077^77'77 


/n3 


I'injinjii 


864 


a* 


/ii277 577^ + 77 (77')^ + 777^77' - 277^77" + 577'' (77')" - 477"77'77" - 3 (77')" 77 


//^2 


//\2 ^ III 


+ ^/i2277 


a 


1677'* + 877 (77')^ + 3277^77' - m^W' + 3477® (77')" + 1877® (77 


/7-1// 


6 - 1 /// 


//\2 


//n3 


-2077" (77')" 77" - 12 {Wf 77" - 1777^77'77 


linj'nj" 


3456 


/22277 77® + 2 {H’f + m^H' + 377" (77')^ 477® - 77"77" + 277 (77')^ - 277'77 


(Al) 


C2 : = 


24 •) 

2 ( 2/2 + 3 / 11 ) (477" + 377') + —/12 1077^ - 37777" + 6 (77') + 1177"77 
^ ^ a"' 

36 

-I- 

,2 


a 

24 


/iii77 (277® - 77") 


+^/22 2077® + 12 (77')^* + 3577^77' - 677®77" - 127777'77" + 3277" (77 
L 

72 . 

-—zfii2H 


t\2 


a* 
144 


877® - 377"77" + 277 (77')^ + 877®77' - m'H 


/i2277 (77" + 277') 1077® - m^H" + 477 {H'f + 477®77' - 677'77 


/22277 (77" + 277')^ 477® - 77"77" + 277 (77')^ - 277'77" 


(A2) 
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C3 := 


-Afio^n + 12/11 - W) + 8/2 - w - 2 uh') 

48 „ 

H— kJi2 


2m^ - m^w + H {u'y + im^w - 2 h'h'' + —/m {2h^ - u^W) 


216 


48 


/22 


60?^' + AH {H'f + im^H' - im'^H 
+22H^ {H'f - 2m^H'H" - 4 {H'f H'‘ 




144 


/iis-H" 2my - IIH^H" + m {H'f + 2m^H' - im'H 


/3-17' 


f^/i22'H" A2W + 20?^ {H'f + Qm^H' - im'^H" + 32H^ {H 


/3 


A0H:^H'H"-28 {H'fH" 


576 


/ 22277 ' 


+ 8{H'f + IAH^H' AH^ - H^H" + 2H {H'f - 2H'H" , (A3) 


//\2 


C 4 : 


= 2 ( 2/2 + 3 / 11 ) (3772 - 77') + —/12 2977^ - 87777" - 5 {H'f + 1077^77' 


36 


—fniH (277^-77") 


r /22 


/2 


mw - 8 (77')'' + 7177^77' - 2077^77" - 287777'77" + 1877' (77 


24 r o 

/ii277 2 877^ - 1177'77" + 677 (77')^ + 2077^77' - 
L 

48 r 

-^/i2277 4277^ + 2077 {H'f + 6877^77' - 1377^77" 


a 


96 


+3277^ (77')^ - 4077'77'77" - 28 (77')' 77 

577^ + 8 (77')^ + 1477'77'1 [477® - 77'77" + 277 (77')^ - 277'77" . (A4) 


//\2 ^ III 


/22277 


16 16 

C 5 := 4 ( 2/2 + 3 / 11 ) 77 + ^/i 2 ( 277® + 77777') + ^/22 [577® + 1477®77' + 877(77')'] (A5) 
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